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Abstract 

We observed the bright unidentified TeV 7-ray source HESS J1616— 508 with the X-ray Imaging 
Spectrometers onboard the Suzaku satehite. No X-ray counterpart was found to a hmiting flux of 
3.1x10"^^ erg s~^ cm~^ in the 2-10 keV band, which is some 60 times below the 7-ray flux in the 1- 
10 TeV band. This object is bright in TeV 7-rays but very dim in the X-ray band, and thus is one of the 
best examples in the Galaxy of a "dark particle accelerator." We also detected soft thermal emission with 
kT ^0.3-0.6 keV near the location of HESS J1616— 508. This may be due to the dust grain scattering halo 
from the nearby bright supernova remnant RCW 103. 
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1. Introduction 

Since the discovery of cosmic rays, the origin of these 
high energy particles has been a mystery. Recent ob- 
servational studies have revealed that some supernova 
remnants (SNRs) have non-thermal X-ray emission (e.g. 
Koyama et al. 1995; Hwang et al. 2002; Bamba et al. 
2003b; Vink, Laming 2003; Bamba et al. 2005), which 
can be interpreted as synchrotron radiation from electrons 
with energies approaching E ^ 10^^ eV. TeV 7-rays have 
also been detected from some of these non-thermal shell- 
type SNRs. Such very high energy (VHE) gamma-rays 
have been explained by either (1) Inverse- Compton (IC) 
upscattering of cosmic microwave background (or other 
lower frequency) photons by the same high energy elec- 
trons giving rise to the X-ray synchrotron emission or (2) 
the decay of neutral pions that originate in collisions be- 
tween high energy protons and dense interstellar matter 



(e.g. Pannuti et al. 2003; Lazendic et al. 2004; Katagiri 
et al. 2005). Thus the combination of VHE gamma- 
ray and non-thermal X-ray emission is expected to be a 
key observational feature of high energy particle accel- 
erators. A high sensitivity survey of the Galactic plane 
in VHE 7-rays has been conducted by the High Energy 
Stereoscopic System (HESS) team, who discovered four- 
teen new VHE 7-ray sources (Aharonian et al. 2005a; 
Aharonian et al. 2006). HESS J1616-508 (hereafter 
HESS J1616) is one of the brightest of these new sources 
and, furthermore, is spatially extended with an angu- 
lar diameter of ~16' in the HESS data. Three known 
objects, the young hard X-ray pulsar PSR J1617— 5055, 
the SNR RCW 103 (G332.4-0.4), and the SNR Kes 32 
(G322.4+0.1), are in the vicinity of HESS J1616 (- 10'- 
15' away from it), but none provides a convincing identi- 
flcation (Aharonian et al. 2005a; Aharonian et al. 2006). 
Before the HESS Galactic plane survey, a few extended 
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7-ray objects with no clear counterpart in other wave- 
bands had been found. These objects have been de- 
noted "dark particle accelerators" (Ubertini et al. 2005). 
They are of interest since they may be indicating the 
acceleration of nucleons; on the other hand, if the TeV 
emission comes from electrons, peculiar conditions such 
as extremely low magnetic fields may be necessary, due 
to the short radiative lifetimes of high energy elec- 
trons (Yamazaki et al. 2006). TeV 2032+410, discovered 
by the High Energy Gamma-Ray Astronomy (HEGRA) 
collaboration (Aharonian et al. 2002), was the first object 
of this kind. In spite of several multiwavelength follow-up 
studies, no compelling counterpart has been found from 
radio to X-ray wavelengths (Mukherjee et al. 2003; Butt 
et al. 2003; Butt et al. 2006). HESS J1303-631 is an- 
other example (Aharonian et al. 2004a). Mukherjee et 
al. (2003) observed this target with the Chandra X- 
ray Observatory for 5 ks, and placed an upper limit of 
<5.4x 10~^^erg s~^ cm~^ to its diffuse hard X-ray fiux in 
the 2-10 keV band. 

HESS J1616 is brighter than both TeV 2032+410 and 
HESS J1303-631 in the TeV regime (see also table 3) 
and is therefore a more suitable case for examining if 
dark accelerators really exist and for finding out just 
how "dark" they are. Since many VHE 7-ray objects 
are X-ray emitters, the X-ray band is the natural place 
to search for a counterpart. To date there have been 
no X-ray observations explicitly targeting HESS J 16 16, 
so we observed HESS J 16 16 with the Suzaku satel- 
lite (Mitsuda et al. 2006). The Suzaku X-ray Imaging 
Spectrometers (XIS) (Koyamaet al. 2006), with their high 
sensitivity and stable low background (Yamaguchi et al. 
2006), are optimal devices to search for diffuse X-ray emis- 
sion. This is especially the case in the hard X-ray band 
(£^>6 keV), where non-thermal emission dominates over 
thermal and the effects of X-ray absorption are minimized. 
By detecting or even just setting a sensitive upper limit 
on the X-ray emission from HESS J1616, it should be pos- 
sible to address the question of what types of particles, 
electrons or protons, are accelerated in the object, and to 
constrain the acceleration mechanism by comparing the 
X-ray and 7-ray emission. In addition to the Suzaku data 
we also used XMM-Newton archival data that partly cov- 
ered the HESS J1616 region. In this work, uncertainties 
are quoted at the 90% confidence level unless otherwise 
stated. 

2. Observations and Data Reduction 

The center of HESS J1616 was observed on 2005 
September 19 during the Science Working Group (SWG) 
phase program. We also observed two nearby positions 
with the same Galactic latitude (BGDl and BGD2 in fig- 
ure 1). These positions were chosen so that no known 
bright X-ray sources fell in the field of view and were 
close enough to HESS J1616 that the brightness of the 
Galactic ridge X-ray emission (e.g. Worrall et al. 1982; 
Koyama et al.. 1986; Yamauchi, Koyama 1993; Sugizaki 
et al. 2001; Revnivtsev et al. 2006) could be estimated 
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Fig. 1. Suzaku field of view overlaid on the HESS smoothed 
excess map in Galactic coordinates (axis labels are in de- 
grees). The circle shows the position and approximate size 
of RCW 103, which is the nearest X-ray bright SNR. The 
scale bar below the figure shows the excess. 

accurately at the source position. The observations are 
summarized in table 1. 

The observations were made with four CCD cameras 
(designated as XIS [0-3]; Koyama et al. 2006) at the foci 
of four X-ray telescopes (XRT; Serlemitsos et al. 2006), 
and a hard X-ray detector (HXD; Kokubun et al. 2006; 
Takahashi et al. 2006). Since the hard X-ray source 
PSR J1617-5055 (Torii et al. 1998) is in the HXD field of 
view, we concentrate on the XIS data. 

The XIS was operated in normal clocking mode during 
all three observations. The edit mode was 3 x 3 or 5 x 5 for 
the HESS J1616 field and we combined the data of both 
modes for our analysis, while only the 3x3 mode was used 
in the background pointings. We used the HE AD AS soft- 
ware version 6.0.4 and version 0.7 of the processed data 
^. Data affected by the South Atlantic Anomaly, Earth 
occultation, and telemetry saturation were excluded. Hot 
and fiickering pixels were also removed. We did not screen 
the data with elevation angle from the bright Earth in or- 
der to maximize the statistics in the hard X-ray band; 
more than 99% of the contamination from the bright 
Earth consists of the O Ka and N Ka lines originating 
in the atmosphere, and the data above 0.6 keV are not 
affected by the contamination. In the following analysis 
we use the XIS data from 0.6 keV to 12 keV. After data 
screening, the effective exposures were 45 ks, 21 ks, and 
24 ks for the HESS J1616, BGDl, and BGD2 regions, re- 
spectively. 

^ Version 0.7 processing is an internal processing applied to the 
Suzaku data obtained during the SWG phase, for the pur- 
pose of establishing the detector calibration as quickly as pos- 
sible. Official data distributed to guest observers are pro- 
cessed Version 1.0 or higher. See Mitsuda et al. (2006) and 
Fujimoto et al. (2006) for more details. 
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Table 1. Log of Suzaku observations. 



Name 


Position (degree) 
/ b 


Start Time 


End Time 




Effective Exposure (ks) 
XISO XISl XIS2 XIS3 


HESS J1616 
BGDl 
BGD2 


332.400 -0.150 
332.000 -0.150 
332.700 -0.150 


2005/9/19 12:09:58 
2005/9/18 22:56:14 
2005/9/20 19:37:50 


2005/9/20 19:37 
2005/9/19 10:59 
2005/9/21 07:28 


50 
02 
46 


45.0 45.4 45.0 45.0 

21.1 23.7 21.2 21.7 
23.5 23.7 23.5 23.5 



3. Analysis and Results 

3.1. XIS image 

We concentrate on the FI CCD data in this imag- 
ing analysis, since non- X-ray background (NXB) dom- 
inates the data of the BI CCD (XISl) in the high en- 
ergy band, especially above 8 keV. In Figure 2, the XIS 
images of the HESS J1616 region are shown for the soft 
(0.6-3.0 keV) and hard (3.0-12.0 keV) energy bands. We 
excluded the corners of the chips that are illuminated 
by the ^^Fe calibration sources, and images from the 
three FI CCDs (XISO, XIS2, and XIS3) were summed. 
Images of the NXB were constructed from night Earth 
data provided by the XIS team and subtracted from the 
HESS J1616 images. After the NXB subtraction, vi- 
gnetting corrections were done as follows. We produced 
simulated XIS images assuming a uniform surface bright- 
ness, using the XRT+XIS simulator xissim ver. 2006-05- 
28. We assumed monochromatic X-rays of energy 1.49 
or 8.05 keV, where the vignetting function is best cali- 
brated (Serlemitsos et al. 2006). The simulated images 
were normalized so that pixels at the optical axis had val- 
ues of ~1. We divided the soft and hard HESS J1616 im- 
ages by the simulated 1.49 and 8.05 keV images. Finally 
the resulted images were rebinned by a factor of 8 and 
smoothed using a Gaussian function with a sigma of 0.^42. 

Since the lower left corner of the field of view marginally 
overlaps the bright SNR RCW 103, corresponding bright 
emission can be seen in the soft band image. On the 
other hand, there is no such structure in the hard band 
image. Furthermore, there is no apparent X-ray structure 
suggesting an X-ray counterpart of HESS J1616 both in 
the soft and hard band images. To examine this more 
quantitatively, we made photon count profiles along the 
strips A A' and BB' with a width of 4' 2 in figure 2, and 
the profiles are shown in figure 3. 

In the soft energy band, the photon count gradually 
decreases along both the X and Y axes, which shows the 
influence of RCW 103. A similar trend can be seen even in 
the hard energy band along the X axis, but we do not see 
it along the Y axis. The trend in the hard energy band 
may have a different origin from the soft band profile: 
the hard X-ray pulsar PSR J1617-5055 or fluctuations of 
the galactic ridge emission may explain the hard band 
profile. We see no systematic trend in either the hard or 
soft X-ray profiles consistent with the TeV 7-ray profile 
of HESS J1616, which can be described by a Gaussian 
function with (j=S'2 (Aharonian et al. 2006). In summary, 
according to our imaging analysis, there is no apparent X- 
ray counterpart to HESS J1616. 



3.2. XIS spectrum 

We made the XIS spectrum of the HESS J1616 region 
from each XIS sensor by extracting X-ray events from 
within a 5' radius of the center of the field of view. We 
also tried to extract spectra for the BGDl and BGD2 re- 
gions in the same way as for the HESS J1616 region, but 
we discovered several X-ray sources in the extraction re- 
gions; there is an X-ray object in the BGDl region (des- 
ignated Suzaku J1614-5114) at (/,6) = (331.°98,-0.°22), 
and two X-ray sources are found in the BGD2 region at 
(/,6) = (332.°67,-0.°19) and (332.°77,-0.°20) (designated 
Suzaku J1617-5044 and Suzaku J1618-5040, respectively). 
According to the SIMBAD Astronomical Database op- 
erated at CDS, Strasbourg, France^, Suzaku J1614-5114 
and Suzaku J1618-5040 are positionally coincident with 
the B-type star CD-50 10270 and the infrared source 
IRAS 16145-5033, respectively. We were unable to iden- 
tify the counterpart to Suzaku J1617-5044. We excluded 
these sources from the BGDl and BGD2 spectra using 2' 
radius circular regions. The spectra from the FI CCDs 
were combined after extraction. 

These spectra contain the NXB. For the most ac- 
curate NXB estimate, we filtered the night Earth 
data so that the cut-off rigidity distribution was the 
same for the HESS J1616, BGDl and BGD2 spec- 
tra (Koyama et al. 2006), and extracted the NXB 
spectra using the same regions in detector coordinates 
(DETX/Y). The NXB spectra thus made were subtracted 
from the HESS J1616, BGDl, and BGD2 spectra, and the 
resulting spectra are plotted in figure 4, where the BGDl 
and BGD2 spectra are renormalized to compensate for the 
difference in the area of the spectral extraction regions. 

Figure 5 is the ratio of the spectrum of the HESS J1616 
region to those of the background regions. The ratio is 
larger than unity below 4 keV, while it is close to unity in 
the 4-8 keV band. The ratio becomes noisy above 8 keV, 
where the NXB dominates the high energy band. 

The most likely origin for the soft band excess in the 
HESS J1616 region is the nearby SNR RCW 103. Here 
we consider whether instrumental effects, i.e., the tail 
of the point spread function (PSF) and stray light of 
the XRT, might be the cause. The HESS J1616 region 
spectrum exhibits emission lines at 0.83 keV, 0.92 keV, 
1.02 keV, and 1.35 keV, which are attributed to Fe XVII, 
Ne IX, Ne X, and Mg XI, respectively, and strongly 
indicate thermal plasma emission. This type of soft 
spectrum resembles that of RCW 103. Using xissim 
we simulated how many photons would fall into the 

^ http://simbad.u-strasbg.fr/Simbad/ 
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Fig. 2. Suzaku XIS images of the HESS J1616 field in the soft (left; 0.6-3 keV) and hard (right; 3-12 keV) energy bands. Lines 
of constant Galactic latitude and longitude are plotted. The scale bars under the figures show the photon count. The images were 
smoothed using a Gauss function with a sigma of 0.^42. Vignetting correction was applied after subtracting non X-ray backgrounds, 
as described in the text. The dotted lines in the left figure show the regions used for the photon-count profiles shown in figure 3. 
The dotted circle with a radius of 5' in the right figure shows the region where we extracted the XIS spectra. 
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Fig. 3. Photon count profiles of the XIS images along the strips AA' (a) and BB' (b) shown in figure 2 The errors are estimated 
at the la confidence level. The curve in the figures shows the Gauss function of cr = 8.^2, which expresses the TeV 7-ray profile of 
HESS J1616. 
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HESS J1616 region from RCW 103, which has a ra- 
dius of and is centered 13' from the HESS J1616 
region. In the simulation, we employed a nonequilib- 
rium ionization plasma model with temperature kT = 
0.3 keV, column density A^h = 7x10^^ cm~^, and ion- 
ization parameter nt=6xl0^ cm~^ yr for the spectrum 
of RCW103 (Gotthelf et al. 1997). As for the metal 
abundances, we set them to 0.5 times the cosmic val- 
ues (Anders, Grevesse 1989). The flux of RCW103 
was set to l.Sx 10~^^erg s~^ cm~^ in the 0.6-2.0 keV 
band (Tuohy, Garmire 1980) and the Chandra ACIS im- 
age of RCW103 (from ObsID 123) was used as the input 
surface brightness distribution. The simulation predicted 
count rates of 1.6x10"^ c s"^ (FI) and 2.4x10"^ c s"^ 
(BI) in the 0.6-2.0 keV band. The measured FI count 
rate from the HESS J1616 spectrum is 4.6x10"^ c s"^ 
with an error of ±0.1x10"^ c s~^ depending on the back- 
ground data set used. The BI count rate is 9.4x 10~^ c s~^ 
with the same uncertainty. These are more than a fac- 
tor of 10 higher than predicted from the instrumental ef- 
fects alone, suggesting that the observed soft emission in 
the HESS J1616 region may have an astrophysical origin. 
However, the current calibrations of the PSF tail and stray 
light at large off- axis angles have large systematic errors. 
A more detailed understanding of the XRT will be needed 
before drawing more definitive conclusions about the pu- 
tative soft emission, and that is beyond the scope of this 
paper. 

Since the spectral ratio becomes noisy above 8 keV, we 
examined the ratio of HESS J1616 to BGDs from 4 to 
8 keV. The average of the ratio and its 99% confidence 
range is found to be l.lOt^ ii and LOltoH in the cases 
of HESS J1616/BGD1 and HESS J1616/BGD2, respec- 
tively. Thus we cannot conclude a positive detection of 
hard X-ray emission from HESS J1616 at the 99% confi- 
dence level with the current statistics of the Suzaku XIS 
data. 

Next we proceed to spectral fitting. Both the FI and 
BI spectra of the HESS J1616 region in the 0.6-8 keV 
band, from which either BGDl or BGD2 were used for the 
background subtraction, were fitted simultaneously. To 
describe the putative soft X-ray emission, we used a ther- 
mal plasma model (the APEC model; Smith et al. 2001), 
and we added a power-law model to quantify any hard 
X-ray emission from the HESS J1616 region. Both com- 
ponents were modified by the same interstellar absorp- 
tion. In the spectral fitting, we used standard RMFs and 
we made ARFs for a flat sky distribution with the soft- 
ware xissimarf gen ver. 2006-05-28. Quoted values for 
fluxes and normalizations in the following are normalized 
to correspond to a uniform sky distribution over a 5' ra- 
dius. The best-fit parameters are summarized in table 2, 
and an example of the best-fitted spectra are shown in 
figure 6. Note that in the spectral fits, the power-law 
component is always statistically significant even at the 
99% confidence level, although it is mainly required to 
fit the spectra in the 2-4 keV band. This energy range 
includes the soft excess discussed above, which may be re- 
lated to RCW 103. If we replace the power-law model by 



Table 2. Best-fit parameters for the model fitting. 



Background^ 



BGDl 



BGD2 



A^H (lO^^cm-^) 
kT (keV) 

Abundance (cosmic) 
/,APE^(0.6-2keV)^ 

normalization (10~^)^ 
xVd.o.f. 



0.50 (fixed) 
7.05x10-^^ 
2.0(fixed) 

7 74+2-73 
' -'^-2.83 

273.11/232 



4 lQ+0.4b 

0.50 (fixed) 
7.55x10 
2.0 (fixed) 

9.00t2"92 

276.45/232 
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a: Data used for background in the spectral fitting. 
b: Flux (erg s~^ cm~^) of the APEC model in the 0.6 - 
2 keV band (uncorrected for absorption), 
c: Photon index of the power-law model. 
d: Normalization of the power-law model defined 
as the differential photon number flux (i.e., pho- 
tons keV~^ s~^ cm~^) at 1 keV. Errors are quoted at 
the 99% confidence level. 

a second thermal component, we obtain fits with similar 
goodness of fit. For example, in the case of the BGD2- 
subtracted spectra, temperatures of kT =0.62 keV and 
kT=2.0 keV plasma with the same metallicity of 0.5 times 
cosmic and Nn =3.5x10^^ cm~^ can fit the spectra with 
X^/d.o.f =265. 50/ 237. The 2.0 keV thermal component, 
which may represent contamination from RCW 103 or be 
part of the putative soft emission, emits 82% of its energy 
flux in the 0.6-2 keV band. Thus the significant detec- 
tion of a hard power-law component does not necessarily 
demonstrate the existence of hard X-ray emission from 
HESS J1616. Given the uncertainties associated with the 
soft diffuse excess we use the power-law component as 
a conservative upper limit to the hard X-ray emission. 
Furthermore the BGD2-subtr acted case gives the larger 
value, and so we take this as the upper limit on the nor- 
malization of the power-law component. Our upper limit 
at the 99% confidence level to the unabsorbed flux in the 
2-10 keV band of HESS J1616 is 3.1xl0-^^erg s-^ cm-^. 

3.3. XMM-Newton analysis 

XMM-Newton (Aschenbach et al. 2000) pointed at the 
X-ray pulsar, PSR J1617— 5055, on 2001 September 3 
(ObsID: 0113050701), and this observation partly cov- 
ered the HESS J1616 region with the EPIC instrument 
(Striider et al. 2001; Turner et al. 2001). We also analyzed 
the XMM-Newton data to check for consistency with the 
Suzaku results. 

The pn camera was operated in timing mode during this 
observation, and the data have no imaging information. 
On the other hand, the MOSl and M0S2 cameras were 
operated in standard full-frame mode using the medium 
filter. We therefore analyzed only the MOS data. We 
used the Standard Analysis System (SAS) software ver- 
sion 6.0.5 for event selection. From the detected MOS 
events, we selected those with PATTERN keywords be- 
tween and 12 as valid X-ray events. Time intervals of 
high and flaring background were rejected by removing 
times when the 10-12 keV count rate was higher than 
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Fig. 4. XIS spectra from the HESS J1616, BGDl and BGD2 regions: (a) combined FI spectra (XIS0+XIS2+XIS3), and (b) BI 
spectra (XISl). Non X-ray backgrounds were subtracted as described in the text. Error bars on the data points are plotted at the 
la confidence level. 



Dark Accelerator HESS J1616-508 



7 



HESS J1616/BGD1 
^ • FI A BI 



(0 



o 



C\2 - 



1 



HESS J1616/:pGD2 
• FI A BI 

4i 




10 



Energy (keV) 



Fig. 5. Ratio of the spectra shown in figure 4. Errors on the 
data points are la. 



> 



o 

I 



: 1 1 1 1 1 

^CTr^E, 







Energy (keV) 

Fig. 6. HESS J1616 spectra of the FI and BI CCDs with the 
best-fit model. The BGD2 data were used for background 
subtraction. 

0.15 c in the full field of view. The resultant exposure 
time was 13 ks for each MOS camera. 

Figure 7 shows the combined MOSl and M0S2 image in 
the 2.0-7.0 keV band. We see no significant excess emis- 
sion at the position of HESS J1616. Because of its higher 
spatial resolution, the EPIC observation is more sensitive 
to point sources than the XIS observation even given the 
shorter exposure time. So we searched for X-ray point 
sources in the combined M0S1+M0S2 2.0-7.0 keV im- 
age using emldetect in the SAS software package. In the 
emldetect software, source significance was determined 
by comparing the numbers of photons in 68% encircled 
energy regions to those in the background map, which 
was prepared using the esplinemap software. The like- 
lihood limit was taken to be 10, corresponding roughly 
to a 4<j detection. The detected point sources, excluding 
those within RCW 103, are shown with plus marks in fig- 
ure 7. Only one point source was detected at the edge 
of the XIS spectral region for HESS J1616, but its fiux 
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Fig. 7. MOSl+2 2.0-7.0 keV image smoothed with a 
Gaussian kernel of cr = 10" . Lines of constant Galactic lat- 
itude and longitude are plotted and labeled in the interior 
of the figure. Plotted are the FOVs of the XIS observations 
(dashed-line squares), the XIS spectral region (thick circle), 
detected point sources (plus marks), and the source and back- 
ground spectral regions (solid and dashed thin circles, respec- 
tively). Projection profile in figure 8 is made along the strip 
AA'. 

in the 2-10 keV band is 2.2xl0~^^erg s~^ cm~^, which is 
smaller than the XIS upper limit for the hard X-ray emis- 
sion from HESS J 16 16 by more than a factor of 10. A 
source was also found at around the center of the BGD2 
region, but we had already excluded it from the XIS anal- 
ysis as Suzaku J1617-5044. 

The detection limit Sum in the 2-10 keV band 
of this search for point sources at the 4a confi- 
dence level was calculated using the formula Sum < 
(4rA/7r6)/(0.68/T), where b is the surface brightness of 
the background map, r is the radius of the 68% encir- 
cled region, / is a counts-to-fiux conversion factor (EOF) 
of 3.8x10^^ cts cm-2 erg-\ and T is the "effective" 
exposure-time corrected for the vignetting. For point 
sources, we obtained Sum =2xl0~^^ erg s~^ cm~^ at the 
center of the XMM field of view, while at the center of 
the XIS field of view, which is away from the optical 
axis. Slim = 3x10""*^^ erg cm~^. 

We made photon count profiles along the strip AA' in 
figure 7 in the 0.6-3 keV and 3-7 keV band (figure 8). 
Note that NXB subtraction and vignetting corrections 
were not applied to the profile. As a reference, we also plot 
the PSF^ at the position of the RCW 103 rim (-6' away 
from the optical axis) plus constant profile in the figure. 
The difference between the soft and hard band profiles 
show an extended soft X-ray halo surrounding RCW 103, 
qualitatively supporting the soft excess seen in the Suzaku 
data. 

We also searched for larger scale X-ray emission at the 

3 In XMM-SOC-CAL-TN-0022 at 

http: / / xmm.esac.esa.int / external / xmm_sw_cal / calib / documentation / index. si 
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Fig. 8. MOSl+2 projection profiles in the 0.6-3 keV and 
3-7 keV bands along the strip A A' in figure 7. 

position of HESS J1616 in the EPIC data. Source and 
background spectra are extracted from the circular regions 
(5' radii) shown in figure 7 (the solid and dashed thin 
circles, respectively). The background region was chosen 
to be at the same galactic latitude as the source region so 
that differences in galactic diffuse emission (Kaneda et al. 
1997) would be minimized. The point sources detected in 
the source and background region were excluded with 2" 
radius circles. MOSl and M0S2 spectra were co- added. 
To avoid the influence of the putative soft emission, we 
ignored the data below 3 keV and fitted the spectrum in 
the 3-7 keV band with a power-law model of a photon 
index F = 2. The 99% confidence limit on the power-law 
normalization is <2. 66x10""^. The XMM- Newton data 
thus set an upper limit of 6.9x10"^^ erg s~^ cm~^ on the 
flux of the hard X-ray emission from HESS J1616 in the 
2-10 keV band. This upper limit is three times larger 
than the Suzaku result. 

4. Discussion 

Neither the Suzaku nor the XMM-Newton data pro- 
vide evidence for hard X-ray emission from HESS J 16 16; 
the XIS sets the most stringent upper limit of 
3.1xl0"^^erg s"^ cm-^ in the 2-10 keV band. The dif- 
ferential photon flux in the TeV 7-ray regime can be rep- 
resented by a power-law of photon index F = 2.35 with a 
total photon flux of (43.3 ± 2.0) x 10"^^ cm-^ s"^ above 
200 GeV (Aharonian et al. 2006). This corresponds to 
an energy flux of 1.7xl0~^^erg s~^ cm~^ in the 1-10 TeV 
band. The ratio of the TeV 7-ray flux to the X-ray flux 
(/xev/Zx) is therefore more than ~55. Table 3 is the 
list of spatially extended VHE objects with X-ray ob- 
servations that we extracted from the H.E.S.S. source 
catalog ^. We also added TeV J2032-F4130 (Aharonian 
et al. 2002) to the hst. Although recent Swift obser- 
vations found possible counterparts of HESS 1614—518 
and HESS J1834-087 (Landi et al. 2006), we do not in- 

4 littp://www.mpi-hd.mpg.de/hfm/HESS/public/HESS_catalog.ht] 



elude them in table 3, since their X-ray data are not 
statistically enough for spectroscopy. Note that the X- 
ray couterparts in the list are not secure. For exam- 
ple, AX J1838. 0—0655 is the most promising candidate 
counterpart of HESS J1837— 069 (Aharonian et al. 2005a; 
Aharonian et al. 2006), it is slightly outside of the HESS 
source extension (Landi et al. 2006). HESS J1616 has the 
most stringent X-ray flux upper limit among the VHE ob- 
jects; yet of more import is its extremely large flux ratio, 
/xev/Zx, compared to the others. Evidently HESS J 16 16 
is a very peculiar object. The bright TeV 7-ray emis- 
sion strongly suggests the presence of some particle accel- 
eration processes there; however HESS J1616 emits lit- 
tle X-ray (or lower frequency radiation) emission and 
thus remained undiscovered prior to the HESS survey. 
HESS J1616 is thus a strong candidate for being a "dark 
particle accelerator." 

Assuming the origin of the TeV 7-ray to be IC scatter- 
ing of the cosmic microwave background by accelerated 
electrons, we can calculate the synchrotron emission from 
these electrons assuming magnetic field values of 5 =10, 1, 
and 0.1 /iGauss (figure 9). The Suzaku flux limit requires 
that the magnetic field in HESS J1616 be less than a few 
micro Gauss, which is the typical interstellar value and far 
below the estimated magnetic field values in other HESS 
sources such as RX J0852-4622 and RX J1713. 7-3946. 
Thus other mechanisms for the TeV emission may be re- 
quired. Aharonian et al. (2006) proposed an asymmetric 
undetected pulsar wind nebula (PWN) powered by the 
young pulsar PSR J 16 17— 505, which is ~13' away from 
HESS J1616, to explain the TeV emission. The spin-down 
luminosity of PSR J1617-505 is 1.6xlO^Vg s"^ (Torii 
et al. 1998), and in this case the luminosity of the PWN 
is expected to be ~10^^ erg s~^ (Cheng et al. 2004) in 
the 2-10 keV band. Assuming the pulsar distance to be 
3.3 kpc (Torii et al. 1998), we obtain an unabsorbed flux of 
7.6xl0~^^ erg s~^ cm~^ in the 2-10 keV band. To recon- 
cile this high flux and the XIS upper limit, we would need 
an absorption of A^h > 2x 10^^ cm~^ for a photon index of 
2. Since the column density is more than 10 times the ab- 
sorption to the Galactic center Baganoff et al. (2003), this 
seems an unlikely scenario to explain X and TeV emission. 
The high /xev/Zx ratio might be explained if HESS J1616 
is an old SNR colliding with a giant molecular cloud, 
since TeV and X-ray emissions are dominated by hadronic 
processes and synchrotron radiation from secondary elec- 
trons, respectively (Yamazaki et al. 2006). However there 
are no radio data indicating the presence of a giant molec- 
ular cloud in the region. In any case, it remains highly de- 
sirable to obtain detailed data of the HESS J1616 region 
in other wavelengths, especially the radio band, to detect 
the synchrotron emission from the accelerated electrons. 

We found extended soft X-ray emission suggesting a 
thermal plasma with /cT~ 0.5 keV both in the Suzaku and 
XMM spectra. Both the temperature and the projection 
profiles of the soft band images (figures 3 and 8) imply a 
close relationship between the soft emission and ROW 103. 
The best-fit TVh of 2-4x10^^ cm-^ in table 2 also sup- 
^port the connection, since these values are close to those 
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Table 3. Spatially extended VHE objects with X-ray observations. 



Name 


Possible counterpart 


TvDe'' 


^ TeV 


PC 

J TeV 




^ X 


fl 


t'~v^\r / t'sr 
J lev / J Js. 


Reference^ 


HESS J0852-463 


RX J0852— 4622 


SNR 


2.1 


6.9 


4 


2.6 


^10 


^ 0.7 


12 3 


HESS T1303— 631 




? 


2.4 


1.0 


20 


2.0 


<0.64 


> 1.6 


4 5 


HESS J1514-591 


PSR B1509-58 


PWN 


2.3 


1.6 


8.6 


2.0 


3.2 


0.5 


6, 7 


HESS J1632-478 


AX J1631. 9-4752 


HMXB? 


2.1 


1.7 


210 


1.6 


1.7 


1.0 


8, 9 


HESS J1640-465 


G338.3-0.0 


SNR 


2.4 


0.71 


96 


3.0 


0.30 


2.4 


8,'l0 


HESS J1713-397 


RX J1713.7-3946 


SNR 


2.2 


3.5 


8 


2.4 


54 


0.065 


11, 12 


HESS J1804-216 


Suzaku J1804-2142 


? 


2.7 


0.48 


2 


-0.3 


0.025 


19 


8, 13 


HESS J1804-216 


Suzaku J1804-2140 


? 


2.7 


0.48 


110 


1.7 


0.043 


11 


8, 13 


HESS J1813-178 


AX J1813-178 


? 


2.1 


0.89 


110 


1.8 


0.70 


1.3 


8, 14 


HESS J1837-069 


AX J1838.0-0655 


? 


2.3 


1.4 


40 


0.8 


1.3 


1.1 


8, 15 


TeV J2032+4130 




? 


1.9 


0.20 


? 


? 


<0.20 


>1.0 


16 


HESS J1616-508 




? 


2.4 


1.7 


4.1 


2.0 


<0.031 


>55 


This work 



a: SNR=supernova reninant, PWN=pulsar wind nebula, HMXB=liigli mass X-ray binary 
b: Photon index of TeV spectra. 

c: Unabsorbed flux in the 1-10 TeV band (in 10~^^ erg s~^ cm~^). 
d: Column density for X-ray spectra (in 10^^ cm~^). 
e: Photon index of X-ray spectra. 

/: Unabsorbed flux in the 2-10 keV band (in 10~^^ erg s~^ cm~^). 

g: (1) Aharonian et al. (2005c), (2) Slane et al. (2001), (3) Our analysis of the ASCA archival data, (4) Aharonian 
et al. (2004a), (5) Mukherjee, Halpern (2005), (6) Aharonian et al. (2005b), (7) DeLaney et al. (2006), (8) Aharonian 
et al. (2006), (9) Rodriguez et al. (2003), (10) Sugizaki et al. (2001), (11) Aharonian et al. (2004b), (12) Slane et 
al. (1999), (13) Bamba et al. (2006), (14) Brogan et al. (2005), (15) Bamba et al. (2003a), and (16) Aharonian et 
al. (2002) 

obtained with the ASCA spectra (Gotthelf et al. 1997). 
The XMM-Newton profile (figure 8) suggests it extends 
up to 15' from the central point source 1E161348-5055, 
which is three times as large as the radius of RCW 103. 
Although we cannot fully reject the possibility, at least for 
the Suzaku data, that this soft emission can be explained 
by an instrumental effect, the presence of this emission in 
the XMM-Newton data argue strongly for an astrophysi- 
cal origin. 

Many bright X-ray sources in the Galactic plane have 
spatially extended X-ray halos, caused by the scattering of 
X-rays into our line of sight by interstellar dust between us 
and the object (Predehl & Schmitt 1995). The fact that 
the hard band profile of XMM-Newton does not extend 
beyond the rim of RCW 103 may support the dust halo 
hypothesis, since the scattering cross section (adust) de- 
pends strongly on photon energy as ddust oc E~'^ (Predehl 
& Schmitt 1995). We will address the soft emission again 
in future when the calibrations have become much robust. 
With careful work it may be possible to model the soft ex- 
tended contamination over the region of HESS J 16 16 and 
further reduce the hard band X-ray flux limit. 

5. Summary 

We observed the bright TeV 7-ray object, HESS J1616, 
with the Suzaku XIS for 45 ks. There is no positive 
detection of hard X-ray emission and we set an upper 
limit of 3.1xl0~^^erg s~^ cm~^ to the 2-10 keV band 
flux. We also analyzed the XMM-Newton data of a 
13 ks observation and obtained a three-times higher up- 
per limit of 6.9x 10~^^erg s~^ cm~^. The deeper expo- 
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Fig. 9. Spectral energy distribution of HESS J1616 from the 
X-ray to TeV 7-ray band. The synchrotron radiation from 
accelerated electrons, which boost the 3 K background up to 
the TeV energy range, is plotted toward the left for three 
different values of the magnetic field. 



sure and low background performance of Suzaku explain 
this difference. The unusually high value of /xev/Zx > 
55 makes HESS J1616 a very peculiar object. We also 
found diffuse X-rays consistent with thermal emission at 
kT ^0.3-0.6 keV extending over this region from the SNR 
RCW 103. This emission is likely the result of scattering 
of X-rays from RCW 103 by interstellar dust. 

The authors are grateful to Profs. W. Hoffman and 
S. Funk for kindly providing us the HESS image. We 
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